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and finally if the density were no greater than that of water the satellite, if 
fluid, would be unstable, the limit being actually within the surface. Since, 
therefore, the satellite could never have existed at its present distance in a 
nebular state, it must follow, if any form of the nebular hypothesis is to be 
accepted, that its original distance was much greater than the present. Can 
we find a probable cause for this ancient disturbance? 

If we suppose the former period of Mars to have been very nearly one- 
sixth that of Jupiter the close commensurability would render the orbit of 
Mars more and more eccentric. The planet in perihelion would thus pass 
through the sun's atmosphere, or rather through the outermost equatorial 
zone of the solar nebula. This resisting medium would not only accelerate 
the motion of Mars but also in a much greater degree that of his extremely 
small satellite. The solar mass contracting more rapidly than the orbit of 
Mars would finally leave the latter moving in an excentric path without 
sensible resistance. 

Other Secondary Systems. — For the first satellite of Jupiter the limit is 
5250 miles, or 4f times the radius of the satellite. For Mimas, the inner- 
most satellite of Saturn, it is less than twice the radius. The rings of 
Saturn, in all probability, could not exist as three satellites, the limits of sta- 
bility being interior to the surface.* 

The effect of perturbation in the dismemberment of comets is known to 
all astronomers. The nucleus of the great comet of 1880, which approach- 
ed within less than 100,000 miles of the sun's surface, must have had a 
density greater than that of granite, as well as a strong cohesive force betw'n 
its parts, in order to withstand the tendency to disintegration during its pe- 
rihelion passage. Had the nucleus been either liquid or gaseous, or even a 
cluster of solid meteorites, the difference between the sun's attraction on the 
central and the superficial parts would have pulled the comet asunder, 
spreading out the fragments into somewhat different orbits, like the meteoric 
streams of August and November. 



LAW OF FACILITY OF ERRORS IN TWO DIMENSIONS. 



BY E. L. DE FOREST, WATERTOWN, CONN. 

In a recent paper (Analyst, Nov., 1880) I gave some account of what 
is known as Lagrange's theorem in probability, and showed how it can be 

*It has been recently shown that Bessel's mass of the ring is much greater than the true 
value. 
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applied to the determination, in the most general manner, of the law of fa- 
cility of errors in one dimension, such as errors in the measured leDgth of a 
line for example. The elementary errors to which any measurement is lia- 
ble, give rise by their combinations to all the errors possible in the sum or 
the mean of several measurements. The sum of k observations is always h 
times their arithmetical mean, so that the sum and the mean follow the 
same law of facility of deviation from their most probable values. The 
probability of a deviation x in the sum is equal to the probability of a devi- 
ation x-i-Jc in the mean. Denoting the quadratic mean error of a single ob- 
servation by r, it was noticed (p. 172) that the quadratic mean error of the 
sum of k observations is r\/h. This is true whether k be finite or infinite, 
being a consequence of my theorem respecting the radius of gyration of the 
coefficients in the expansion of a polynomial to the k power. (Analyst, 
May, 1880, p. 75.) Hence the quadratic mean error of the arithmetical 
mean, for any value of h, is strictly equal to 

ri/h-T-h = r-^^/k. 

It will be borne in mind that r here is in strictness the square root of the 
mean of the squares of the deviations of all the possible values of a single 
observation, from their arithmetical mean, and that in forming this arith- 
metical mean, and the mean of the squares of the deviations, each possible 
value, and its squared deviation, is taken with a weight proportional to the 
probability of its occurrence. 

If the number of observations is large, the determination of the law of 
facility of error reduces itself to finding the limiting form of the series of 
coefficients in the expansion of a polynomial to an infinitely high power. 
This, as I had before shown, is the probability curve 

y - -yz e ■ (!) 

We will now consider the somewhat analogous case of errors in two di- 
mensions, or errors in the position of a point in a plane. Taking the true po- 
sition of the point as an origin, and drawing through it any rect. axes, we 
refer to these axes the positions of the other points corresponding to the sev- 
eral elementary errors to which the observations are subject. Let Ax and 
Jy denote the units of measure in the x and y directions. If one elementary 
error acting aloDe would place the observed point at the position x = 2Jx, 
y = 3dy, for example, we will say that k 2t 3 is the probability of the occur- 
rence of that error ; and so on for the other errors. Now write the proba- 
bilities X as coefficients, and their sub-indices as exponents of f and r], in the 
polynomial, 

1 £— m n m I 2 E m n m 
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For convenience, only the four corner terms are written here, as in a pre- 
vious article of mine (Vol. VII, p. 41). The omitted terms can be readily 
constructed from these, since all sub-indices, exponents and coeffic'ts which 
vary from term to term, are understood to vary by differences equal to imi- 
ty. The greatest elementary error, positive or negative, which can occur, 
is supposed not to exceed m units in either the x or y directions. The poly- 
nomial, if written in full, would form a square, with 2m -f- 1 rows and 2m 
-fl columns, containing (2m+l) 2 terms in all. 

If any particular elementary error cannot occur, its probability X is of 
course zero ; and any polynomial can be put in the square form by adding 
terms, if necessary, with zero coefficients. By an extension of the principle 
of Lagrange's theorem, which a little consideration will render obvious, it 
appears that when k observations are made, the probability that the sum of 
the errors in the x direction will be tAx, and that at the same time the sum 
of those in the y direction will be vAy, is the coefficient of ? 'rf in the expan- 
sion of the polynomial (2) to the k power. The same coefficient is also the 
probability that the mean of the a>errors will be t-t-k units and the mean of 
the ^-errors v~-k units. The units of measure Ax and Ay can be taken as 
small as we please, and the exponents in the polynomial will be regarded 
as whole numbers. 

In k observations, the most probable total effect, in the x and y directions 
of the elementary error 2 Ax and 3 Ay for example, will be 
kX^ z (2Ax), kX^ z (ZAy), 

and so on for other errors. The most probable algebraic sums of the x and 
y errors then will be 

kAx(- mX - m '™ \ mL - m ), kAy( ml -™>™ L™k-_ \ ,<n 

\—mX_ m< _ m | mX m< _V \— mX_ mt _ m | — mA m , _J ' w 

so that these are approximately the exponents of $ and y in that term in the 
expansion of (2) whose coefficient is a maximum. Hence the most probable 
error in the mean of the k observations, in the x and y directions respect- 
ively, is found by dividing the expressions (3) by k. It will be seen that 
the quotients thus obtained are the statical moments of the coefficients X 
about the axes of Y and X, these coefficients being regarded as masses in a 
system of material points. The same quotients also represent the lever arms 
of the system about those axes, since the sum of all the masses X is unity. 
From the principle established in my article of March '80, p. 46, it follows 
that the lever arms of the system of coefficients in the expansion of (2) to the 
k power, about the axes of X and Y respectively, are k times what they are 
in the first power. According to (3) therefore, the maximum coefficient in 
the expansion will coincide approximately with the centre of forces, or cen- 
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tare of gravity, of the whole expanded series of coefficients, whether k be fi- 
nite or infinite. We shall hereafter get in another way a proof that this is 
true at the limit, when k becomes infinite. 

Now to find the probability of a given error in the mean result of h 
observations, we must be able to find, approximately at least, the coefficient 
of any term of given position in the expansion of the polynomial (2) to the 
k power. To this end, we shall proceed to demonstrate a relation which 
exists between the coefficients in any square group of (2m+l) 2 terms in the 
expansion, and the co-ordinates which fix the position of the group. 
This is equivalent to finding a "multinomial formula" for two variables, 
analogous to the known formula for one variable which was employed by 
me in Analyst, Sept., 1879, p. 141. That formula is demonstrated most 
easily by the method of indeterminate coefficients. (See for instance Ber- 
trand's Calcul Diffe'rentiel, p. 330.) A similar process will be followed here, 
and we shall then pass to the differential equations of the surface which is 
the limiting form of the coefficients in the expansion of (2), in a manner 
somewhat analogous to that by which, in my own article just cited, the dif- 
ferential equation of the limiting curve was reached. 

The coefficients in the expansion of a polynomial such as (2) to the k 
power are not altered by adding a constant number to all the exponents of 
£, or to all those of rj, or to both together. If then we write 

u = «Q.n? 7 n KngT (4) 

«o.o£VKo£V' {) 

and make n = 2m, the whole number of terms will be the same as in (2), 
and if we also make all the coefficients a equal to the coefficients X which 
are in the corresponding positions, the coefficients in the expansion of (4) to 
the k power will be the same as those in the expansion of (2). Denoting 
them by B, the expanded polynomial may be written 

the number of terms in it being (kn+lf. We now have 

u k = U. (6) 

Differentiation with respect to ? and tj gives 

*»'-'© = % *»*-($ - f • < 7 > 

and dividing these by (6) and clearing of fractions, we get 

»"(&-•(£)• w ®-(w)- (8) 

By differentiation of (4) and (5) with respect to £, we find 



du la^J'Tf \mt n J"-^vj n dU li} 1Jt ,fV' 1 knB^ J P-iif m 
d$ la 1 , o eVI«a*.o**~ 1 7 ' dS 15 ll0 fV | *»»"»*». 5*— 1 7° 

and by substitution in the first of the two equations (8), 

ft/ -Bo,*.C l*l-B*.,fe.C*V' \ / l«i..CV h«n,„C"~V" \ * ~) 

"" Voo.oCV I a^oCVAl-Bi. oCV I *»**.. oC*"" 1 ," J J 
Also differentiating (4) and (5) with respect to //, and substituting in the 
second of the equations (8), we get a second equation analogous to (10), 
which I omit, to save space. The reader can easily construct it for himself. 
Each member of (10) contains the product of two rectangular polynomials. 
If we multiply together any two such polynomials, for instance 

/ ct 0> xv k.c ;"" \/ oo..cy Kw:y° \ nu 

WoCV I ^oCV/Wocy l °..oC>v' { } 

their product will also be a rectangular polynomial, and any term in it will 
have for its coefficient the sum of products of a and c, such that this sum 
can be expressed in the rectangular form. For example, the coefficient of 
£ p fi q in the product of (11) is 

a 0,q c p,0 I a P.q°0,0 w /io\ 

**o, °p. q I ®p, o c o, a 
Now forming in this way the coefficient of C -1 /*' m * ne product in each 
member of (10), and equating the two to each other by the principle of in- 
determinate coefficients, we get 

j/ ra rt0 B 0t ,\ la^oB^yA = ra 0i ,B r , \ la (r _^ s B li0 

W-,.-B ,oi l"'l,«-#(r-l), o' ra 0.oA, «,) l«(r-l). 5 llS ' ' 

Likewise forming the coefficients of CM' -1 in the products in the omitted 
equation similar to (10), and equating them to each other, we have 

^ A a f,l-"0,(»-1) I l a 0,l-Br,(t-l) \ _ l a 0.(»-l)-ffr,l | lq,,(,-i)Z? .i /, .. 

Vsa rjJ j5 0i0 \sa 0tS B rt0 I sa 0>0 2? riS jsa Pi0 B 0t ,' ^ ' 
These two equations express a relation between the coefficients a and B 
in the rectangular groups 

a o, > ; a r, > •£><>, s : B r< , . 

'„"'■'„ > R R '' \ lt} ) 

u 0, : u r, JJ 0, ; JJ r, 

the one group containing coefficients in the first power of the given poly- 
nomial, while the other group contains an equal number of coefficients in 
its expansion to the k power. Since r and s are arbitrary, suppose each of 
them to be greater than n. This will not alter the value of the given poly- 
nomial, nor that of its expansion, provided we suppose the coefficients a to 
be equal to zero whenever either of the sub-indices of a is greater than n. 
*For want of sorts J and fi are here written to represent Xi and Eta, respectively. 
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Omitting therefore in (13) and (14) all terms which contain such values of 
a as factors, we shall have remaining 

J, / "«», Q-Bfr-n),, | 1«1 , Q-B(r- 1 ),« \ ra(l,n-Br, (,_„) | (?*-»)«„, n-^(r-,.),(.-n ) ' 

Uo.^^,.,,)^..,,) | la j ,„£(,_!), (,_„)/ ra 0>0 B ri , j (r-^)a Bt0 £ (r _) iS 

]c( "•! -"(T-tt), (»-l) | 1 Q 0, 1 -" r, (»-l) \ = _ l CT 0,n-"r, (»— n ) | (^~ w )Q». n-P(r-n),(»— n) 

^"n, i ^M, (m) | ""|, . -Or, (m) ' Sa oa B r< , | SCl n>0 -D(r_n), « > 

(J= s _n). (16) 
These two equations express a relation between the coefficients a and i? in 
the square groups 

a<),n ' a n> „ -0(r— ii),« ' -"r, < /I 7\ 

,----, p- -.- p , (XI) 

a 0, : a »,0 -"(r-n),^— n) ; -"r, (i-«) 

the first group containing the (n + 1 ) 2 coefficients of the given polynomial 
(4), while the second contains an equal number of coefficients in its expan- 
sion to the k power, these last being situated so as to form any square group, 
such that the highest sub-index of B in the x direction is r, and the high- 
est in the y direction is s. The equations (16) can be advantageously mod- 
ified as follows; in the first members add a column or a row, with coeffic'ts 
zero ; and reverse the positions of the terms in the second members, as if by 
turning the two rectangular tables round through 180° in their own plane. 
We thus get 

k I na,,, o.B ( ,-,,), . 1 0a , „B r< , \ _ (r-n)q„ i0 g (r _„ ) . , [ (r-0)a B r , , ^ g 
V na n> B -B (r _ B) , ( ,_ B) | 0a 0i n B r< ( _ n) / (r-n)a„ t B -B (r _ B) , ( ,_ B) | (r-0)a 0i „J3,. j( ,_ B) ' 
while the second equation, being formed in a similar way, may be omitted 
here. In the second member of (18), let that portion which does not have 
the coefficient r be transferred to the first member, and we get 

(k-j-i)( na "' <f ~ n) ' * I ° ' ° r - - ^ = r[ a ". «('•-'')• ; 1 a o.o-"'-, ; ) n 9) 

\na B „X>( r _ n )(,_ B ) | Oa 0-B i> rj (j _ B ) / \tf„ t n-D(r-n),(»-fi) I a 0,n-"r, (»-n) ' 

the second equation being omitted as before. Since the coefficients in the 
polynomial (2) are represented by those in (4), and we have n=2m,the num- 
ber of terms in each row or column is 2m -f- 1 in the first power of (4), and 
2km + 1 in the expansion to the h power, so that one coefficient in the ex- 
pansion will be in the middle of the square. This construction does not 
affect the generality of our proof of the limiting form of the expansion of a 
polynomial, for any given polynomial can be made to take this square form 
with one term in the centre, by adding terms on one or more sides, if need 
be, with coefficients zero. Now let us fix the position of the square group 
(17) of coefficients B, not by the extreme ranks r and s reckoned from the 
place of the first coefficient JB , but by the ranks i andj of the middle term 
of the group, reckoned from the place of the middle coefficient in the whole 
expansion, in the x and y directions respectively. This gives 



r = m{k+l)-+i, * = m(k+l)+j (20) 

Substituting these values for r and s in (19) and its companion equation, 
and transposing from the second members to the first ones those portions 
which have the coefficient m(£+l), the former equation becomes 

(A-l-lH mO "< < r -">. " ~ ma 0,q£>r, s ^ t '( q «.0-"(i--n) | g 0,0-"r,« \ /OJS 

\WM n , M i>( r _ „),(s-k) — ma o,n-&r, (e— u) ' ^ a n, n-0(r— >i),(s-n) | a 0,nB r ,(s— n)/ 

and I omit the second one as before. Let us now change the notation so that 
instead of the coefFs a and B as in (17), we shall write X and £, as follows, 

^ — ffi, m ' "m,m ^(i—m),(j-{-m) \ ^(i+m),(j-fm) fOON 

r , r \aa) 

"—m,—m ; "m,—m v (i-m).(j—m) ; v (i-\~m) t (y— m) 

Substituting these for the a and B in (21) and its companion equation, and 
using V as an auxiliary letter, we get the following results, 

^m, ~my(i— m),0'+m) | ^~m,—m^(i+m) ,(J-{-m) | 

— 1 7 r~j 1 ' 

"m, m ' J (i—m),(3'—m') \ ^—m t m (i+tfi),(/ — m) 

-m'j i +m),(J+m ) __ / * \ y i (23) 



Wjm, ~m'(i—m),(j+m) | m ^-m,-mt( . +m),(/+pi ) __ / * \ xr ! 

m X mm f(j_ m ) j y_ m ) I Wl/_ OT . im t(i+m),(j— m) \A-f-l/ | 



^m, — mt(»— m).Q"+in) | WM— m,-M^(i+m),CH-» i) ( 7 j TT 
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These two equations, like (16), express the general relation between the 
(2m+l) 2 coefficients in any polynomial of square form, and an equal num- 
ber of coefficients forming a square group in any part of the expansion of 
the polynomial to the k power. It amounts to a "multinomial formula" 
for two variables. Its analogy to the formula for one variable, as given by 
me in ANALYSt, Sept., 1879, p. 142, is apparent. The position of the sq're 
group of coefficients I in the expansion is fixed by the ranks i and j of its 
middle term reckoned from the middle coefficient of the whole expansion, 
so that, referring the polynomial and its expansion to the same axes, the co- 
ordinates of the middle coefficient of the group will be 

x = %Ax, y=jAy- (24) 

It will be seen that the terms in (23) are made up of products of X into I, 
and the particular ones to be multiplied together are readily found by rota- 
ting the X table in (22) through 180°, in its own plane, and then applying 
it to the I table, and multiplying together the terms which coincide. The 
sum of all the products thus formed is V. And it will be seen that if these 
products are regarded as parallel forces acting perpendicularly to the plane 
of XT, and the intervals between them are taken as unity of distance, their 
moments about vertical and horizontal axes through the middle of the group, 
constitute the first members of the two final equations in (23). The rela- 
tion here proved is a general one, holding true for any given polynomial, 
whether its coefficients X are positive or negative. 

[To be continued.] 



